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Abstract 

The Plebahski-Demiahski metric, and those that can be obtained from it by tak- 
ing coordinate transformations in certain limits, include the complete family of 
space-times of type D with an aligned electromagnetic field and a possibly non-zero 
cosmological constant. Starting with a new form of the line element which is better 
suited both for physical interpretation and for identifying different subfamilies, we 
review this entire family of solutions. Our metric for the expanding case explicitly 
includes two parameters which represent the acceleration of the sources and the 
twist of the repeated principal null congruences, the twist being directly related to 
both the angular velocity of the sources and their NUT-like properties. The non- 
expanding type D solutions are also identified. All special cases are derived in a 
simple and transparent way. 
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1 Introduction 

The complete family of type D space-times with an aligned non-null electromagnetic field 
and a possibly non-zero cosmological constant A can be represented by a metric that 
was given by Plebahski and Demiahski £Q together with those that can be derived from 
it by certain transformations and limiting procedures. These solutions are characterised 
by two related quartic functions whose coefficients are determined by seven arbitrary 
parameters which include A and both electric and magnetic charges. For the vacuum 
case with vanishing cosmological constant, they include all the particular solutions that 
were identified by Kinnersley j2J. For the sub-cases in which the repeated principal null 
congruences are expanding, these metrics have been analysed further by Weir and Kerr 
where the relations between the different forms of the line element can be deduced. 
(They have also been given independently by Debever and Kamran 0] and Ishikawa and 
Miyashita 0.) 

Unfortunately, many particular and well-known type D space-times are not included 
explicitly in the original form of the line element. They can only be obtained from it by 
using certain degenerate transformations. Moreover, the parameters that were introduced 
in the original papers are not the most useful ones for a physical interpretation of the 
solutions. One purpose of this paper is therefore to present a new form of the Plebahski- 
Demiahski metric in which the parameters are given a clear physical meaning and from 
which the various special cases can be obtained in a more satisfactory way. In this way 
we will clarify the complete family of solutions. 

A new look at these solutions seems to us to be particularly important at the present 
time because they are now being used in new ways. People working in semi-classical quan- 
tum gravity have used these metrics to investigate the pair production of black holes in 
cosmological backgrounds (see e.g. [H] [E])- Others are working to extend these solutions 
to higher dimensions. Yet these solutions are still not well-understood at the classical 
level of general relativity. In particular, the physical significance of the parameters em- 
ployed in the original forms are only properly identified in the most simplified special 
cases. Thus the emphasis here is not to look for the most general metric form which 
covers all cases. This was achieved many years ago with the work of Carter j^j, Debever 
et al |10j . Garcia D. jTTj and others. Rather, our purpose is to cast the metric in a form 
in which the parameters employed have clear physical interpretations and thus to classify 
the complete family of solutions in a way that clarifies their physical properties. 

In order to achieve this, we will first modify the original Plebahski-Demiahski metric to 
include two parameters, denoted by a and a>, which respectively represent the acceleration 
of the sources and the twist of the repeated principal null congruences. Then we will 
derive a more general form of the line element which explicitly contains all the well-known 
special subfamilies, at least for the cases in which the repeated principal null congruence 
is diverging and the orbits of the Killing vectors are non-null. (The case in which the 
group orbits are null is not considered here.) 

The solutions contained in this family are characterized by two generally quartic func- 
tions whose coefficients are related to the physical parameters of the space-time. One 
of our main purposes is to clarify the physical meaning of these coefficients which have 
traditionally been misinterpreted in the general case. In particular, we locate the cos- 
mological constant in the most appropriate place and identify the relation between the 
Plebahski-Demiahski parameter n and the NUT parameter. It has also been traditional 
to use an available coordinate freedom to remove the linear term in one of these quartics. 
However, as noted by Hong and Teo in a particular case ^2] > [Sj , it is much more helpful 
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to use the freedom to simplify the roots of this quartic. It will be shown that this helps 
very significantly both in the interpretation of these solutions and also in performing the 
associated calculations. 

In section 3, we will review the family of solutions of this type which possess space- 
like surfaces of positive curvature. These include the solutions which are considered to 
have greatest physical significance as they relate to the fields around generalised forms of 
black holes. After that, we consider the complete families of non-accelerating and then 
non-twisting solutions in more detail. 

In practice, it is convenient to distinguish the cases in which the repeated principal 
null directions are either expanding or non-expanding. A separate general form of the 
metric is required for the non-expanding case. Such a form is derived in section [7| using 
a degenerate transformation of our initial metric. As the twist also vanishes in this limit, 
it follows that these solutions must be type D members of the solutions of Kundt's class 
- a result that is demonstrated explicitly. 

Finally, we identify the members of this family of solutions which reduce to the 
Bertotti-Robinson solution, or other direct product space-times, in an appropriate limit. 



2 The Plebahski— Demiahski metric 

We consider here a general family of type D solutions of Einstein's equations including a 
generally non-zero cosmological constant A. These may be vacuum or include a non-null 
electromagnetic field such that the two repeated principal null congruences of the Weyl 
tensor are aligned 1 with the two principal null congruences of the non-null electromagnetic 
field. 



2.1 The original form of the metric 

Let us start with the Plebahski-Demiahski metric PQ (see §21.1.2 of [15] ) which is given 
by 



ds 2 



where 



Q(dr - p 2 da) 2 V(dr + f 2 da) 2 f 2 +p 2 , ^ f 2 +p 2 

; ; ; ; dp dr 

X— pr) 2 y r 2 + p 2 r 2 + p 2 V Q 



V = k + 2np - ep 2 + 2mp 3 - (k + e 2 + g 2 + A/3)p 4 , 
Q = (k + e 2 + g 2 ) - 2mf + er 2 - 2nf 3 - (k + A/3)f 4 



(2) 



and m, h, e, g, e, k and A are arbitrary real parameters 2 . It is often assumed that rh and 
h are the mass and NUT parameters respectively, although this is not generally the case. 
The parameters e and g represent electric and magnetic charges. The parameter 7 that 
is used in [T] and [12] is obtained by putting k = j — g 2 — A/6. However, as shown in 
[T7] , it is more convenient for physical interpretation to include the cosmological constant 
in the form given in (J2J) . 

1 Other non-null type D electrovacuum solutions exist in which only one of the principal null congru- 
ences of the electromagnetic field is aligned with a repeated principal null congruences of the Weyl tensor 
(see for example |14|V However, such solutions are not considered here. 

2 Debever [16^ had previously found a type D metric with an equivalent set of seven arbitrary param- 
eters. 
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As is required for type D space-times of this type, the general family of solutions 
represented by ((TJ) admits (at least) two commuting Killing vectors d„ and d? whose 
orbits are spacelike in regions with Q > and timelike when Q < 0. Type D solutions 
also exist 3 in which the group orbits are null, but these are not considered here. 

2.2 A modified form of the metric 

For purposes of interpreting the Plebahski-Demiahski metric, it is convenient (see 
and |21]) to introduce the rescaling 



P 



'atop, 



a 



3 a > 



T 



(3) 



with the relabelling of parameters 



m + in 



(-)' 

y 0J J 



[m + in) 



01 , V 

e + ig = -{e + ig), 

OJ 



e=-e, 

OJ 



k = a 2 k. (4) 



This introduces two additional parameters a and oj. With these changes, the metric 
becomes 



ds 2 



1 



apry 



Q 



r' z + oj 2 p 2 



(dr — up 2 da)' 



P 



r' z + oj 2 p 2 



(udr + r 2 daf 



r 2 + oj 2 p 2 l2 r 2 + oj 2 p 2 2 
dp dr 



P 



Q 



(5) 



where 



P = P(p) = k + loj^np - ep 2 + 2amp 3 - [a 2 (oj 2 k + e 2 + g 2 ) + cu 2 A/3]p 4 
Q = Q(r) = (oj 2 k + e 2 + g 2 ) - 2mr + er 2 - 2aoj~ 1 nr 3 - (a 2 k + A/3)r 4 , 



(6) 



and m, n, e, g, A, e, k, a and oj are arbitrary real parameters of which two can be chosen 
for convenience. It should be emphasised that, apart from A, e and g, the parameters 
included in this metric do not necessarily have their traditional physical interpretation. 
They only acquire their usual specific well-identified meanings in certain special sub-cases. 
Adopting the null tetrad 



rv 



apr) 



2(r 2 + oj 2 p 2 ) 
(1 — apr) 



r 2 d T — ojd c 



Qd r 



2(r 2 
1 - 



F oj 2 p 2 ) 
apr) 



^=(r 2 d T - ojdaj + y/Qd, 



(7) 



2(r 2 + u 2 p 2 ^ 



ojp 2 d T + d a ) +iy/Pd 1 



3 Metrics in which the group orbits are null can be obtained by using a degenerate coordinate trans- 
formation. This is described in detail by Garcia and Plebahski ^Sj (see also p322 of JS]) wn0 confirm 
that a particular case of Leroy is included in this general family. All such type D space-times with 
null group orbits have been given by Garcia and Salazar |2()j . A generalized form of the metric Q which 
includes all known cases of this type has been given by Debever, Kamran and McLenaghan 
although this form is not well suited for the interpretation of the solutions. A form of the metric which 
covers the cases of null and non-null orbits simultaneously was used by Garcia ^Jj. (See also [2*2^.1 
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the spin coefficients are given by 



K=a=X=u=0 
p = H = 



Q (l+iaup 2 ) 



2(r 2 + uj 2 p 2 ) (r + iuip) 

2\ 



P 



[u — tar 



T = IX 



6 = 7 = T 



2(r 2 + uj 2 p 2 ) (r + iuip) 

(1 — apr) 



Q 



4 V 2(r 2 + u 2 p 2 ^ 



Q' 



a = j3 = - 



P 



4 V 2(r 2 + u 2 p 2 ) 



-2ap- (1- apr) — 
(r + tujp) Q 



(1 — apr) P' 

2co- + 2iar + i(l — apr) — 

(r + iup) y ' P 



These indicate that the congruences tangent to Z M and n M are both geodesic and shear- 
free but have non-zero expansion. It can be seen that the twist of both congruences is 
proportional to w. In some particular cases u is directly related to both the angular 
velocity of sources and the effects of the NUT parameter (see 

In terms of the tetrad (J7J), the only non-zero component of the Weyl tensor is 



^2 = —(m + in) 



apr 



r + iuop 



+ (e 2 + g 2 ) 



1 — apr \ 1 + apr 
r + iup ) r — iojp 



(9) 



This confirms that these space-times are of algebraic type D, and that the tetrad vectors 
/ M and n M as chosen above are aligned with the repeated principal null directions of the 
Weyl tensor. The only non-zero components of the Ricci tensor are 



1 



(e 2 + g> 



(1 — apr) 4 " 
(r 2 + u 2 p 2 ) 2 



A 



NP 



I A. 



(10) 



Together, these indicate the presence of a curvature singularity at p = r = which, if 
contained within the space-time, may be considered as the source of the gravitational 
field. 

Having introduced a and u as continuous parameters, we are free to use the rescaling 
(J3J) with (jlj) to scale the parameters e and k to some specific values (without changing their 
signs). For example, we could set them to the values +1, or —1, but it will generally be 
more convenient to scale them to some other appropriate values. It is clear that e and g 
are the electric and magnetic charges of the sources and A is the cosmological constant. 
For certain choices of the other parameters, it is found that m is related to the mass of 
the source and n is related to the NUT parameter (although it should not be identified 
with it in general). 

To retain a Lorentzian signature in (jSJ), it is necessary that P > 0. And, since P(p) 
is generally a quartic function, the coordinate p must be restricted to a particular range 
between appropriate roots. If this range includes p = 0, it is necessary that k > and, 
in this case, it is always possible to use the scaling (jHJ) to set k = 1. However, this is not 
necessary, and we simply note that the requirement that P > may place some restriction 
on the possible signs of the parameters e and k. In fact, various particular cases can be 
identified and classified according to the number and types of the roots of P and the range 
of p that is adopted (see e.g. [Sj). Points at which P = generally correspond to poles 
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of the coordinates. By contract, surfaces on which Q = are horizons through which 
coordinates can be extended. It is also significant that 

Q(r) = -a\*P(±)-U^ + A. 

ar 6 \a z J 

When A = 0, this relates the number and character of the roots of these two quartics. 
However, when A ^ 0, this correspondence is obscured. 

The line element © is flat if m = n = 0, e = g = and A = 0, but the remaining 
parameters e, k, a and uj may be non-zero in this flat limit. Moreover, it is not immediately 
obvious that the metric © includes the Schwarzschild-de Sitter solution, the Reissner- 
Nordstrom solution, the Kerr metric, the NUT solution or the C-metric, which are all 
known to be of algebraic type D, or the Robinson-Trautman type D space-times. In the 
following sections, it will be shown that a simple transformations of (j3J) leads to a form 
which explicitly includes all these well-known special cases. 



3 A general metric for expanding solutions 



When A = 0, the line element (jSJ) already contains the Kerr-Newman solution for a 
charged rotating black hole. It also contains the charged C-metric for accelerating black 
holes. However, it does not include the type D non-singular NUT solution [23]. To cover 
all these cases and their generalizations, it is necessary to introduce a specific shift in 
the coordinate p. In fact, this procedure is essential to obtain the correct metric for 
accelerating and rotating black holes. We therefore start with the metric © with (JBJ), 
and perform the coordinate transformation 



V 



I a _, 
- + -P, 

UJ UJ 



T = t- 



(l + a) 



UJ 



a = — - 



a a 

where a and I are new arbitrary parameters. By this procedure, we obtain the metric 



(11) 



ds 2 



02 \f 2 [dt-(a(l-f) + 2l(l-p) 



P r 



adt - (r 2 + (l + a) 




(12) 



where 



n 

p 2 
P 

Q 



1 (1 + ap)r . 

uo 

■r 2 + (l + apf , 



a + aip + a 2 p 



7$ 



a 3 p 3 



a 4 p*, 



{ui 2 k + e 2 + g 2 ) — 2mr + er z 



2a-r 3 - 

UJ 



a 2 k + ±)r* 
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and we have put 



a i 



a 2 



a 3 




a 

2a— m - 

UJ 

2 

a I 2/ i 2 i 2 

— (u k + e + g 

.UJ A 



a. 



uj and A. 



These solutions generally have seven essential parameters m, n, e, g, 
They also have two parameters k and e which can be scaled to any convenient values. 
In addition, we have the further parameters a and I which can be chosen arbitrarily. In 
practice, it is convenient to choose a and / to satisfy certain conditions which simplify the 
form of the metric, and then to re-express n and uj in terms of these parameters. 

The properties of the solutions in this family depend critically on the character of the 
function P(p). In fact, as an arbitrary quartic, P can have up to four distinct roots, and 
Lorentzian space-times only occur for ranges of p for which P > 0. When more than one 
such range exists, the different possibilities correspond to distinct space-times which have 
different physical interpretations. 

When P has no roots, this function can only be positive and p G (— oo, oo). 

For the cases in which P has at least one root, without loss of generality we can choose 
the parameters a and I so that such a root occurs at p = 1. The metric (|12|) is then regular 
at p = 1 which corresponds to a coordinate pole on an axis, and it is then appropriate to 
take 6 periodic coordinate. 

When another distinct root of P exists, it is always possible to exhaust the freedom in 
a and I to set the second root at p — — 1. The metric component a(l — p 2 ) is then regular 
at this second pole while the component 21(1 — p) is not. Thus, the metric is regular 
at p — 1, but a singularity of some kind occurs at p — —1. (In fact, unless I = 0, the 
region near p — — 1 contains closed timelike lines.) With this choice, and for the positive 
curvature case in which both poles are located on a continuous axis, it will be shown 
that a corresponds to a Kerr-like rotation parameter for which the corresponding metric 
components are regular on the entire axis, while I corresponds to a NUT parameter for 
which the corresponding components are only regular on the half-axis p = 1. 

We have now introduced through (fTT|) a shift and scaling of p such that, if P has at 
least two roots, then it admits the two factors (1 — p) and (1 + p). Thus 

P= (1 - p 2 )(a - a 3 p - a 4 p 2 ), 

which implies that the above coefficients must satisfy the conditions 



a i + a 3 — 0, Oq + Q-2 + &4 = 0. 



(13) 



These conditions provide two linear equations which specify the two parameters e and n 
in terms of a and I as 

■• 2 A' 



— + 4a- m - a 2 + 3Z 2 ) 

a 1 — l l uj 



a 



to 



(oj 2 k + e 2 + g 2 ) + 



(14) 



n 



UJ 2 kl 

a 2 -I 2 



UJ 



a " ' ,2,„ , 2 , 2- A 



(uj 2 k + e 2 + g 2 ) + 



(15) 
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Equation (JT5|l explicitly relates the Plebahski-Demiahski parameter n to the NUT pa- 
rameter I. With these definitions, we then obtain that 

a = ^ k L - 2a- m + 3a 2 ^(uj 2 k + e 2 + g 2 ) + I 2 A. 
a 2 — I 2 u uj 2 

The character of the solution then partly depends on whether this expression for ao 
is positive, negative or zero. If it is non-zero, the scaling freedom can then be used to 
set ao = ±1. This equation then effectively defines the parameter k. Thus, for any given 
value of ao, the constant k is given by 




The original metric (JJ) contained the three parameters h, e and k. In the above 
argument, we have increased the number of such parameters to include n, e, k, a, uj, a 
and I, but we have also introduced three constraints that are effectively represented by 
(|14j) . (|15|) and One remaining (scaling) freedom is therefore still available which may 
be used to set uj to any convenient value (assuming a and I do not both vanish). There 
thus remain the three parameters a, a and I in addition to m, e, g and A. 



4 Generalised black holes (P has two roots, ao = 1) 

In this section, we will concentrate on the physically most relevant particular case of the 
line element (|12|) for which P has at least two distinct roots and ao > 0, so that we can 
set a = 1. In this case, the surfaces spanned by p and have positive curvature. 



4.1 The complete family of black hole-like space-times 



For the case considered here, p is taken to cover the range between the roots p = ±1 and 
it is natural to put p = cos 9, where 9 G [0, tt]. In this case, the metric (|T2j) becomes 



where 



and 



ds 2 



J_(Q 

fi 2 \p 2 



dt - {a sin 2 9 + 4/ sin 2 



Pi 

Q 



dr 2 



P 



adt 



(r 2 + (a + l) 2 )d<jj 



tsin 2 9d9 2 



n 

P 2 

p 

Q 



a 



1 (I + a cos 9) r 

UJ 

r 2 + (I + acos9) 2 

sin 2 9 (1 — a 3 cos 9 — a 4 cos 2 9) 



{uj 2 k + e 2 + g 2 



a 3 



a 

2a— m 

UJ 

, a 2 



2mr + er 
al 



n o 
2a-r 3 

UJ 



4a z 



UJ 



(uj z k + e z + g 1 



2, A 1 

a k + — r 



A 
A-al 



a 4 = —a —(ujk + e + g 



A 



;i7) 



18) 



(19) 



uj* ' 3 
with e, n and k given by ()14j ) - (fT^|) . It is also assumed that | a-3 1 and \a^\ are sufficiently 
small that P has no additional roots with 9 G [0, tt]. This solution contains eight arbitrary 
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parameters m, e, g, a, I, a, A and u>. Of these, the first seven can be varied independently, 
and io can be set to any convenient value if a or I are not both zero. 

It was shown in [21J that, when A = 0, the metric (|17jl represents an accelerating and 
rotating charged black hole with a generally non-zero NUT parameter. However, an arbi- 
trary cosmological constant is now included so that the background is either Minkowski, 
de Sitter or anti-de Sitter space-time. For the vacuum case in which A = and e = g = 0, 
the general structure of this family of solutions is given in figure ^ The special cases are 
generally well-known and will not be discussed here. 



mlaa (a>l) 



Accelerating and Rotating 
Black Hole with nut 



mlaa (l>a) 



Accelerating NUT 
solution with rotation 







a = ov 


a = 


maa 




mla (a>l) 


^ mla (l>a) 


Accelerating and 
Rotating Black Hole 




Kerr solution 

with NUT 


NUT solution 
with rotation 



a = 




ma > 


< 


C-metric 



a = o 




Figure 1: The structure of the family of solutions represented by ((T7I) when A = 0, e = g = and 
m/0. This family has four parameters m, I, a and a. An accelerating Kerr solution with a small 
NUT parameter has been distinguished from an accelerating NUT solution with a small rotation 
as their singularity structures differ significantly even though their metric forms are identical. 
For the same reason, the Kerr-NUT solution has similarly been divided. An accelerating NUT 
solution without rotation has not been identified. All the special cases indicated have obvious 
charged versions and versions with a non-zero cosmological constant. 



The non-zero components of the curvature tensor are given by (jHJ) and (fTTH) in which cop 
is replaced by I + acos#. It can be seen that the metric (|T7j) has a curvature singularity 
when p 2 = 0. If |/| < |a|, this occurs when r = and cos# = —l/a. On the other 
hand, if |Z| > \a\, p 2 cannot be zero and the metric is non-singular. These two cases 
have to be considered separately as they clearly have very different global and singularity 
structures. This distinction is also indicated in figure [T] 

When a 2 > I 2 , the metric (fTTjl has a curvature singularity at r = 0, cos# = —l/a. 
And, since the space-time is asymptotically flat at conformal infinity where Q = 0, the 
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range of r must be given by r 6 (0, r^) where 

'"oo 



if a cos 9 > —I 



a(l + acos6) (20) 
oo otherwise 

In the alternative case in which a 2 < I 2 there is no curvature singularity and the range of 
r can also take negative values. However, the coordinate r does not extend to conformal 
infinity in all directions (i.e. for all values of 9) when a^O. When I > \a\, r G (—oo, r^), 
where = ua" 1 ^ + acos9)~ 1 . And, when I < — \a\, r 6 (— | r oo|)°°)- In these cases, the 
space-time is not asymptotically flat at infinite values of r, as this does not correspond to 
conformal infinity. This is a natural feature of accelerating coordinates. 

As fully described in [21], conical singularities generally occur on the axis. However, 
by specifying the range of appropriately, the singularity on one half of the axis can be 
removed. For example, that on 9 = tt is removed by taking e [0, 2tt(1 + a 3 — a^ 1 ). In 
this case, the acceleration of the "source" is achieved by a "string" of deficit angle 

4ira 3 

d o = t- (21) 

1 + CI3 — CI4 

connecting it to infinity. Alternatively, the singularity on 9 = could be removed by 
taking e [0, 2n(l — a 3 — azi) -1 ), and the acceleration would then be achieved by a 
"strut" between the sources in which the excess angle is given by 

An a 3 , . 

-5„ = - 3 — . 22 

1 — a 3 — 04 

The expressions (J21|) or (J22j) are closely related to the tension or stress in the string or strut 
respectively and these should be equal for any given acceleration a, at least according to 
Newtonian theory. However, the deficit /excess angles are the same fractions of the range 
of the periodic coordinate in each case. This, presumably, corresponds to an equality of 
forces in the general case. 

It should also be recalled that, when / ^ 0, the metric (fTTj) has an additional singularity 
on 9 = 7r which corresponds to the "axis" between two causally separated "sources". 
However, this can be switched to the other half-axis by the transformation t' = t — 4/0. It 
can thus be seen that the topological singularity on the axis which causes the acceleration, 
and the singularity on the axis associated with the NUT parameter and the existence of 
closed timelike lines (see e.g. P6j), are mathematically independent. They may each be 
set on whatever parts of the axis may be considered to be most physically significant. 

Intriguingly, there exists a possibility that the conical singularity may vanish on both 
halves of the axis simultaneously. This occurs when 03 = 0: i.e. when 



2a 2 l(u 2 k + e 2 + g 2 )-\u 2 lK 



auom. 



However, this generally 4 only occurs when the NUT parameter / is non-zero and, in this 
case, one of the "half-axes" is singular in a different way and is surrounded by a region 
containing closed timelike lines. Thus, the presence of charges or a cosmological constant 
is not sufficient to cause the sources to accelerate: some string-like structure is required 
(at least within this family of solutions). 

4 It may be observed that this condition for a completely regular axis may also be satisfied when a = 0, 
/ = uj = and 2a(e 2 + g 2 ) = m. This is the special case of the extremely and oppositely charged 
C-metric. 
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Some special cases will now be considered in more detail. Namely, the non-accelerating 
solutions and those with no NUT parameter. The general case with A = has been 
described in detail in [23] for which the charge-free case is summarised in figure ^ Their 
particular subcases are well known. For the case with a = 0, a is a redundant parameter, 
indicating that there are no accelerating NUT solutions without rotation. 



4.2 The Kerr-Newman-NUT-de Sitter solution (a = 0) 

When a = 0, (JTBJ) becomes u 2 k = (1 — / 2 A)(a 2 — I 2 ) and hence (fT3j) and (fTB]) become 

e= l-(§a 2 + 2/ 2 )A, n = I + §(a 2 - 4/ 2 )/A. 

The metric is then given by ()17j) with 



n 



r 2 + (7 + a cos I 



P = sin 2 0(1 + |AaZcos# 



3 



e 2 + S 2 ) 



+ ±Aa 2 cos 2 0) 
2mr + r 2 — A 



/ 2 )/ 2 + (ia 2 + 2/ 2 )r 2 + 



l r 4 
3' 



This is exactly the Kerr-Newman-NUT-de Sitter solution in the form which is regular on 
the half-axis 6 = 0. It represents a non-accelerating black hole with mass m, electric and 
magnetic charges e and g, a rotation parameter a and a NUT parameter I in a de Sitter 
or anti-de Sitter background. It reduces to known forms when I = or a = or A = 0. 

In this Kerr-Newman-NUT-de Sitter solution, it is important to distinguish the two 
cases in which | a \ is greater or less than | / 1 , as indicated in figure ^ When a 2 > I 2 , 
k > 0, the metric has a Kerr-like ring singularity at r = 0, cos 8 = —l/a, and the 
range of r is given by (J2(Jj) . This case represents a Kerr-Newman-de Sitter solution (a 
charged black hole) with a small NUT parameter. Alternatively, when a 2 < I 2 , k < 0, the 
metric is singularity free, and the range of r includes negative values. This case is best 
described as a charged NUT-de Sitter solution with a small Kerr-like rotation. Although 
these two cases have identical metric forms, their singularity and global structures differ 
substantially. 



4.3 Accelerating Kerr— Newman— de Sitter black holes (/ = 0) 

In the physically most significant case in which a is arbitrary but / = 0, (fTp^l implies that 
uj 2 k = c 
hence 



It is then convenient to use the remaining scaling freedom to put uo = a, and 



2/ 2 , 2 , 2\ 

a (a + e + g J 



|Aa 2 



k 



n = —aara. 



For this case, it can be seen explicitly that the Plebahski-Demiahski parameter n is 
non-zero, while the NUT parameter / vanishes. The metric (|T7j) now takes the form 



ds2 = ^ll 



dt — a sin 2 6 1 



o 2 P 
Q P 2 



adt - (r 2 + a 2 )d(p 



tsm 2 6d9 2 



(23) 
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where 

Q = 1 — ar cos 9 



p 2 = r 2 + a 2 cos 2 9 

P = sin 2 6»(l - 2am cos 9 + [a 2 (a 2 + e 2 + g 2 ) + §Aa 2 ] cos 2 5 
g = (( a 2 + e 2 + g 2 ) - 2mr + r 2 ) (1 - a 2 r 2 ) - ±A(a 2 + r 2 )r 2 . 
The only non-zero components of the curvature tensor are given by 




r — ia cos 9 I \ r + ia cos 9 



3 



1 / 2 2x (1 _ ar COS6 1 ) 4 

$n = o e + y rrr 2 — 2^ and A - 

These indicate the presence of a Kerr-like ring singularity at r = 0, 9 = |. Thus, we may 
take the range of r as 



r G (0, Too) where 



or 1 sec if #<tt/2 
oo otherwise 



When A = 0, the metric corresponds precisely to that of Hong and Teo [13] (and 
described in [2JJ which represents an accelerating and rotating black hole without any 
NUT-like behaviour and in which the acceleration is characterized by a. In this case, if 
m 2 > a 2 + e 2 + g 2 , the expression for Q factorises as 



Q = (r_ — r)(r + — r)(l — a r ) 

where 



r± = m ± J m 2 — a 2 — e 2 — g 2 . (24) 

The expressions for r± are identical to those for the locations of the outer (event) and 
inner (Cauchy) horizons of the non-accelerating Kerr-Newman black hole. However, in 
this case, there is another horizon at r = a^ 1 which is already familiar in the context of 
the C-metric as an acceleration horizon. For the case in which A ^ 0, the locations of the 
horizons are modified. 

It may also be observed that 

a 3 = 2am, a 4 = — a 2 (a 2 + e 2 + g 2 ) — |Aa 2 , 

so that the deficit angle of the string causing the black holes to accelerate, or the excess 
angle of the strut between them, is obtained immediately using (}2~T]) or (}2*2"j) respectively. 
There are no closed timelike lines near the axis, confirming that this is the appropriate 
metric for describing a pair of accelerating and rotating black holes 5 as argued by Hong 
and Teo |TS] for the case when A = 0. 

The metric (J23j) nicely represents the singularity and horizon structure of an accelerat- 
ing charged and rotating black hole in a de Sitter or anti-de Sitter background (see |30j). 
It represents the space-time from the singularity through the inner and outer black hole 



5 It should be emphasised that this solution, which has no NUT-like properties, is different from the 
one that is obtained by putting n — and which is usually called the "spinning C-metric". This case, 
whose properties have been described in [57]-|2nj etc. (at least when e = g = and A = 0) still retains 
NUT-like properties such as the existence of closed timelike lines near one half of the axis. 
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horizons and out to and beyond the acceleration horizon. However, it does not cover the 
complete analytic extension inside the black hole horizon. For this, Kruskal-Szekeres-like 
coordinates are required as described for example for particular cases in [HI] . Neither does 
it cover the complete analytic extension beyond the acceleration horizon. This requires 
a transformation to boost-rotation-symmetric coordinates as has been given for the gen- 
eral case with A = in [32]. This shows that the complete space-time actually contains 
two causally separated charged and rotating black holes which accelerate away from each 
other in opposite spatial directions. 



4.3.1 The charged C-metric with a cosmological constant 

For the case in which a = 0, the metric ()23|) reduces to the simple diagonal form 



(1 — otcos6>) 2 \r 2 Q P J 

where 

P = sin 2 e(l- 2am cos 6 + a 2 (e 2 + g 2 ) cos 2 #) , 

Q = ( e 2 + g 2 - 2mr + r 2 )(l - a 2 r 2 ) - |Ar 4 . 

When A = 0, this is exactly equivalent to the form for the charged C-metric that was intro- 
duced recently by Hong and Teo [O] using the coordinates x = — cos 8 and y = — 1/ (ar). 
It describes a pair of black holes of mass m and electric and magnetic charges e and g 
which accelerate towards infinity under the action of forces represented by a conical sin- 
gularity, for which a is precisely the acceleration. In this case, the acceleration horizon 
at r = a -1 is clearly identified. However, when A / 0, the location of all horizons is 
modified. 

Further properties of the charged C-metric in a (anti-)de Sitter background have been 
analysed in [17] and However, these worked with a different form of the metric 

to that above. 



5 Twisting but non-accelerating solutions (a = 0) 

Let us now return to the general family of expanding type D solutions given by the metric 
f!12[) but concentrate here on the particular case for which a = 0, implying that Q = 1. 
The line element is thus 



p 2 



&t - (a(l - p 2 ) + 21(1 - p) 



P 

7^2 



adt 



(r 2 + (/ + a) 2 )d0] 2 -tdp 2 , 



(25) 



where 



r 2 + + ap) 2 



P = a + dip + a 2 p 2 + a 3 p 3 + a 4 p 4 

Q = (u 2 k + e 2 + g 2 ) - 2mr + er 2 - |Ar 4 , 
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and 

a = a - 2 (uj 2 k + 2nl - el 2 - §AZ 4 ) 
a x = 2a- 1 (n -el- §AZ 3 ) 

a 2 = -e - 2M 2 (26) 
a% = — |Aa/ 
a 4 = — |Aa 2 . 

As explained in section El the properties of the solutions in this family depend critically 
on the character of the quartic P(p). However, it may be noticed here that the physi- 
cal parameters m, e and g do not appear explicitly in this function, but the remaining 
parameters n, k, e, A, a, I and u must be such that there exists at least one range of p 
in which P > 0. When more than one such range exists, the different possible ranges 
correspond to distinct space-times which have different physical properties. 

Although it is generally possible to choose / to put ai = 0, it is more convenient to 
choose I to simplify the factors of P. In fact, when Aa = 0, P is a quadratic which can be 
put into one of its distinct canonical forms. Otherwise, P can be expressed as the product 
of two quadratics, one of which can be similarly simplified. 

5.1 Non-accelerating cases with A = 

When both a = and A = 0, P reduces to a quadratic and we obtain the line element 
(J25|) in which 

p 2 = r 2 + (I + ap) 2 
P = a + dip - ep 2 
Q = (uo 2 k + e 2 + g 2 ) — 2mr + er 2 , 

with 

ao = a~ 2 (uj 2 k + 2nl — el 2 ^ , a\ = 2a~ x (n — el) . 

In this metric, a and I can be chosen for convenience and, in addition, we still have 
two scaling freedoms. The resulting family of space-times generally has five arbitrary 
parameters m, n (preferably expressed in terms of /), e, g and a (replacing uS). 

It is now possible to fix a and I (equivalent to the linear transformation ()11|) in p) to 
put P (which is necessarily positive) into one of its six distinct canonical forms, namely 

P: 1-p 2 , p 2 + l, p 2 -l, p 2 , p, I. 

When the electromagnetic field vanishes, this yields the family of class II solutions iden- 
tified by Kinnersley j2J. Let us now consider each of these possibilities explicitly. 

When e > 0, we can put e = 1 and set I = n so that a\ = 0. We can then set 
uj 2 k = a 2 — I 2 , so that ao = 1, and thus P = 1 — p 2 . This case corresponds exactly the 
Kerr-Newman-NUT solution described in subsection 14.21 above, in which it is natural to 
put p = cos 8. In the vacuum case, it is the Kerr-NUT solution, which is the Kinnersley 
class II. A solution. 

When e < 0, we can put e = —1 and set I = —n so that again a\ = 0. We 
can then set u 2 k = I 2 + eoa 2 , where eo = 1,-1,0, so that P = p 2 + eo and 
Q = (e 2 + g 2 + I 2 + eoa 2 ) — 2mr — r 2 . For the vacuum case and with different values 
of e , these correspond to the Kinnersley class II. B, II. C and II. D solutions respectively. 
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(When eo = 1 it is convenient to put p = sinh^, and when eo = — 1 it is convenient 
to put p = cosh 9.) 

Finally, we consider the case in which e = 0. When n ^ 0, it is convenient to set 
a = 2n and u 2 k = —2nl so that a\ = 1 and a = 0. With these choices, the 
metric functions become P = p and Q = (e 2 + g 2 — al) — 2mr. For the vacuum case, 
this is the Kinnersley class II. E solution. Alternatively, for the case in which n = 0, it is 
convenient to set uo 2 k = a 2 so that = 1 and the metric functions become P — 1 
and Q = (e 2 + g 2 + a 2 ) — 2mr. For the vacuum case, this is of Kinnersley class II. F. 

5.2 Non-accelerating cases with A ^ and n = 

When both a = and n = 0, it is appropriate to set I = so that P is a quadratic in p 2 . 
In this case, we obtain the metric (|25|) with 

p 2 = r 2 + a 2 p 2 

P = a~ 2 u 2 k - ep 2 - |Aa 2 p 4 

Q = ( e 2 + g 2 + u) 2 k) - 2mr + er 2 - |Ar 4 . 

Various particular cases now need to be considered. For different values of the essential 
parameters, the free parameters have to be chosen to obtain a convenient range of p in 
which P > 0. Here, we will only consider A ^ 0, since all the cases with A = are already 
described in the previous subsection. 



5.2.1 The Kerr— Newman— de Sitter solution 

When k > we may choose u 2 k = a 2 and then set e = 1 — |Aa 2 provided the signs 
are consistent. Thus we obtain 

P=(l-p 2 )(l + Ua 2 p 2 ) 



Q = ( e 2 + g 2 + a 2 ) - 2mr + (1 - |Aa 2 )r 2 - 
= (e 2 + g 2 ) - 2mr + (r 2 + a 2 )(l - |Ar 2 ) 



|Ar 4 



which is obviously the I = subcase of the space-times discussed in subsection 14.21 The 
metric can then be written in the standard form of the Kerr-Newman-de Sitter solution 
in Boyer-Lindquist-type coordinates [SZj, [SHj using the transformation 



p = cos 9, 



t = tH _1 , 



(27) 



where S = 1 + ^Aa 2 . This puts the metric in the form 



ds 2 = 
where 



A r r - 2 -~l 2 

dt — a sin 6 dcf) 



E 2 p 2 



Af> sin 2 9 r 



E 2 p 2 



adt - (r 2 + a 2 )d0 



2 -^ldr 2 -^ 2 



Aq 



p 2 = r 2 + a 2 cos 2 9 , 



A r = (r 2 + a 2 )(l - \Ar 2 ) - 2mr + (e 2 + g 2 



A f , 



1 



|Aa 2 cos 2 9 . 



(28) 



(29) 



Formally, there is no need to introduce the constant scaling S in t and <fi. However, this is 
included so that the metric has a well-behaved axis at 9 = and 9 = n with e [0, 2tt). 
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It should also be noted that, if A < 0, P(p) = has four roots: p = ±1 and 
p = -. If ^Aa 2 > —1, the above solution is valid for < 9 < it. However, if 

r V — A a 3 ' — — ' 

|Aa 2 < — 1, the solution will only be valid for the smaller range of 9 for which A# > and 
the space-time in this case does not represent a black hole in an anti-de Sitter background. 

5.2.2 Other cases 

When e = and A > 0, it is necessary that k > 0. In this case, p varies between the 
two real roots of P. However, when A < 0, k can take any value and p will take various 
appropriate ranges which all extend to ±oo. 

When e ^ and k > 0, we may choose u 2 k = a 2 . The Kerr-Newman-de Sitter 
solution is obtained as above when the sign of e is such that it is possible to make a 
rescaling to set e = 1 — |Aa 2 . However, if the signs are inconsistent, it is always possible 
to use the scaling to set e = — l + |Aa 2 . In this case, we obtain P = (l+p 2 )(l — |Aa 2 p 2 ). 
If A < 0, this is a valid solution for all p G (— oo, oo). However, if A > 0, it is valid for 

When e ^ and k < 0, we can choose u 2 k = —a 2 and set e = — 1 — |Aa 2 provided 
the signs are consistent. In this case P = (p 2 — 1)(1 — |Aa 2 p 2 ). We can then choose 
p = cosh 9, so that A r = -(r 2 + a 2 )(l + |Ar 2 )-2mr + (e 2 +fi' 2 ) and A e = l-|Aa 2 cosh 2 9. 
However, even this is only valid over the entire range of 9 if A < 0. If < A < 3a" 2 , P 
generally has four roots and two distinct ranges of p are possible. And if A > 3a" 2 , we 
can again put p = cos 9. 

If the signs are inconsistent in the above case for which e ^ and k < 0, we can set 
e = 1 + |Aa 2 . In this case, P = —(1 +p 2 )(l + |Aa 2 p 2 ). Of course, this is only positive 
if A < and p 2 is sufficiently large. 

Finally, we must consider the remaining case in which e 7^ and k = 0. In this case, 
P = —p 2 (e + |Aa 2 p 2 ) which is only permitted if e and A are not both positive. If at least 
one of these parameters is negative appropriate ranges of p can easily be obtained. 

5.3 Non-accelerating cases with A 7^ and a = 

Let us now consider the complementary case in which both a = and a = 0. However, 
the coefficients of P given in (j2B|) must be bounded in the limit as a — * 0. We thus set 

u 2 k + 2nl - el 2 - |A/ 4 = a a 2 
n — el — |AZ 3 = \a\a 
e + 2M 2 = -a 2 

0, so that 03 = = a 4 . Thus, the metric function P in 

P = a + aip + a 2 p 2 . 

e = -a 2 - 2A1 2 , 
n = —a 2 l — f A/ 3 , 
uj 2 k = a 2 l 2 + M\ 



and then take the limit a 
is again the quadratic 

We also obtain 
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and the metric (J23j) can be written as 

ds 2 = Q(dt-2l{l-p)d(f>) -^-{r 2 + l 2 )(Pd<p 2 + ^dp 2 ) (30) 

Q P 

where 

a 2 {l 2 - r 2 ) + (e 2 + g 2 ) - 2mr - A (|r 4 + 2/ 2 r 2 - I 4 ) . 

In this form, and with appropriate values of do that will be clarified below, it may immedi- 
ately be recognised that the three solutions with different values of a 2 correspond exactly 
to the general family of NUT solutions with additional non-zero electric and magnetic 
charges and a non-zero cosmological constant. 

We now consider all possible distinct subcases according to the usual canonical forms 
of P in which the coefficients a, take the values 1, —1 or 0. 




5.3.1 The Taub-NUT-de Sitter solution 

In the case in which ao — 1, ax — and a 2 = — 1 so that P = 1 — p 2 and \p\ < 1, we can 
put p = cos 9 so that the metric (j3*Uj) can be written as 

ds 2 = Q( dt _ Al sin 2 | # ) 2 _i!_ (r 2 + /2) ( dQ 2 + gin 2 q ^ (31) 

Q 

where 

r 2 _ 2mr - I 2 + e 2 + g 2 - A (|r 4 + 2/ 2 r 2 - Z 4 )" . 

This is exactly the known extension of the Taub-NUT metric to include non-zero charges 
and a cosmological constant in which / is the NUT parameter. It is contained in the 
case described in subsection 14.21 for a = 0. 




5.3.2 Other cases 

There are five remaining cases to be considered. However, it is generally convenient here 
to set t = t' + 21(f). 

For the case when a — 1, a± — and a 2 = 0, in which P — 1, the metric (J3Ti|) becomes 

,2 dr 2 



ds 2 = Q(dt' + 2lpd(f) 



?r - l r 2 + Z 2 )(d0 2 + dp 2 ). 



where 



Q 



r 2 + p 



-2mr + e 2 + g 2 - A (|r 4 + 2/ 2 r 2 - V 



(32) 



(33) 



For the case when ao = —1, a\ = and a 2 = 1, in which P = p 2 — 1 and \p\ > 1, we 
put p = cosh i? so that the metric (|30J) becomes 



ds 2 = g(dt' + 2/cosh J Rd0 



2 dr 2 



+ l 2 )(dR 2 + sinh 2 i?« 



where 



l 

r 2 + / 2 



_ r 2 _ 2mr + I 2 + e 2 + g 2 - A (|r 4 + 2/V - / 4 ) 



(34) 



(35) 
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However, in this case, it would have been more convenient to revert to the coordinate t. 
With this, the metric would be regular at R = 0. This would then corresponds to an axis, 
so that (ft may be interpreted as an angular coordinate. 

The above two cases, together with the case ao = 1, a\ = = — 1 that was 
described in the previous subsection, constitute the three members of the general family 
of NUT (Kinnersley's class I) solutions when e = g = A = 0. They thus generalise these 
solutions to include non-zero electric and magnetic charges and a non-zero cosmological 
constant. They correspond, in the order given above, to non-rotating (a = 0) subcases of 
the solutions of Kinnersley's classes II. A, II. F and II. C generalized to include charges and 
a cosmological constant. However, for the family of solutions being discussed here, there 
are apparently still three additional cases to be considered. 

Let us now consider the case in which ao = 0, a\ — 1 and a 2 = 0. In this case, P = p, 
and the metric takes the form 

ds 2 = Q(dt' + 2lpd(f>Y - ^C- - (r 2 + Mpdcf 2 + - dp 2 
v 1 Q v P 

where Q is given by (|33|). Putting p = jp 2 and = 2ip, the metric becomes 

ds 2 = Q(dt' + /p 2 d0) 2 - - (r 2 + / 2 )(dp 2 + p 2 d^ 2 ) . (36) 

Q 

This is a non-rotating (a = 0) solution of Kinnersley's class II. E with additional parame- 
ters. However, the further transformation t' = r + Ixy, p = \/x 2 + y 2 , ip = tan _1 (7//a;) 
transforms ()36|) to a form which is identical to that of (}3*2j) above. Thus, in this non- 
rotating limit, solutions of Kinnersley's classes II. E and II. F are equivalent even with 
additional charges and a cosmological constant. 

For the case when a 2 = 1, ai = and a = 0, in which P = p 2 , we put p = e x so that 
the metric (|3T)|) becomes 

ds 2 = Q(dt' + 2Ze*d0) 2 - - (r 2 + Z 2 )(d X 2 + e 2 *d0 2 ) (37) 

Q 

where Q is given by (|35|). This is a non- rotating (a = 0) solution of Kinnersley's class 
II. D with additional parameters. However, the transformation 

sin tp 

e x = cosh R + sinhi? cosip, <fi 



together with 



cos ip + coth R 
l ( tanh(-|) + cos tp\ 



simp 



t' = T + 2ltfj + 41 tan 

takes the metric (JH7j) exactly to the form (J21J). 

Finally, for the case when a-i — 1, a\ — and a = 1, in which P = 1 +p 2 , we put 
p = sinh ip so that the metric (j3*Uj) becomes 

ds 2 = Q (df + 21 sinh V d<j>) 2 - -Z- - (r 2 + I 2 ) (d^ 2 + cosh 2 i> d0 2 ) (38) 

Q 
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where Q is given by (|35|). This is a non- rotating (a = 0) solution of Kinnersley's class 
II. B with additional parameters. However, applying the transformation 

= sinhx + cosh x cosh as, = r— -—. -. — , 

cosh a + tanh \ 

i, , , , i /l + tanh(|) coshcA 

t' = r - 2/ tan" 1 sinh a) + 41 tan" 1 ^ 

\ smh a J 

to the metric (|37j) takes it exactly to the above form (J38|) . This confirms that these 
generalised solutions of Kinnersley's classes II. B, II. C and II. D (i.e. with a 2 = 1 and any 
value of do) become completely equivalent in the class I solutions for which a = 0. 

It is therefore concluded that the family of charged NUT-de Sitter solutions which 
have five continuous parameters m, /, e, g and A come in just three different varieties 
according to the values of the discrete parameter a 2 = ±1, 0. Of these, the case in which 
a 2 = — 1 may be referred to as the charged Taub-NUT-de Sitter solution. 

Finally, let us remark that when I — 0, e = g — and A = (m ^ 0), the three 
possible solutions are the so-called A- metrics |4U) . When a 2 = —1, it is the AI metric 
which is the Schwarzschild solution. When a 2 = 1, it is the All metric which is its 
hyperbolic equivalent. And when a 2 = 0, it is the AIII metric which is either the type D 
Kasner solution or Taub's plane symmetric solution according to the sign of m. 



6 Accelerating but non-twisting solutions (uj = 0) 

In this section we will consider the twist-free Plebahski-Demiahski solutions ©, © in 
which uj = 0. In this case it is necessary that n = 0, but n' = uj'^-n may be non-zero and 
the metric takes the diagonal form 



1 — apr) 2 




(39) 



where 



P = k + 2n'p — ep 2 + 2amp i — a (e + g )p , 

1m e 2 + a 2 ( 40 ) 
q = e - fH + U-L - 2an'r - (a 2 k + A/3)r 2 , 

and m, n', e, g, A, e, k and a are arbitrary real parameters. (It is particularly significant 
here that the cosmological constant does not appear explicitly in the expression for P.) 
These solutions are in fact the type D Robinson-Trautman solutions. This follows from 
the observation that the repeated principal null congruences are geodesic, shear-free, 
twist-free and expanding. 

When a = 0, P is quadratic and canonical coordinates can be used. For example, 
when e > 0, we can take e = 1, remove n' and set k — 1. It is then natural to put 
p = cos 9 so that P = sin 2 6 and the resulting metric is the familiar form of the 
Reissner-Nordstrom-de Sitter solution. For alternative canonical choices the equivalent 
zero and negative curvature solutions are obtained. These cases in which a = are the 
/ = subcases of the family of solutions described in section 15.31 

When a ^ 0, we note that P is generally a quartic. As such, it can always be written 
as the product of two quadratic factors. However, we can now consider the transformation 
which leaves the conformal factor in (|3*^j) invariant, namely 



(41) 
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where j3 is a constant, Under this freedom, the metric (J39j) is unchanged with e, g and A 
also unchanged, but the other parameters are transformed as 



m = m — 2a/3(e + g ) , 

(42) 



e = e — Qaf3m + 6a 2 /3 2 (e 2 + g 



n > = n'-Pe + 3a(3 2 m - 2a 2 /3 3 (e 2 + g 2 
k = ~k + 2pti - (3 2 e + 2a(3 3 m - a 2 (3 4 {e 2 + g 2 ) . 

Unless e = g = and e 2 < 12am/', it is possible to use the freedom (|41|) to remove the 
parameter n' . However, it is more useful in general to use this freedom, together with a 
rescaling of coordinates, to simplify one of the quadratic factors of P(p). In particular, 
we can introduce a new parameter e' and use the available freedoms to set 

, ka 2 (e 2 + q 2 ) . kam 

£ = e , — , ri = -, 

e e 

so that 

, 2 / 2am a 2 (e 2 + g 2 ^ 



P=(k- e'p 2 ) 1 —p + 



e' 



Q=[l-—r"\[r'-~, | ±r 2 




where e' and k can be scaled to appropriate values. 

Consider first the case in which e' and k are both positive and scaled to unity: i.e. 
when 

e = 1 — a 2 (e 2 + g 2 ), k — 1, n' — —am. 



In this case 



p 2 ) (l - 2amp + a 2 (e 2 + g 2 )p' 



2m i e 2 + g 2 \ A 



Q = (1 - a 2 r 2 ) 1 - — + ^—^ - r 2 



which is precisely the form of the charged C-metric given in [24] except that an arbitrary 
cosmological constant is now included (see also section E.3.1JI . This case thus represents 
a charged black hole which is accelerating in a de Sitter, anti-de Sitter or Minkowski 
background with acceleration a, under the action of a string-like structure represented by 
a conical singularity. 

Other space-times are obtained when k and/or e' are negative, or k = 0, although 
these may have no physically significant interpretation. For these cases, the parameter a 
may not have the interpretation as an acceleration. As above, the range of p must always 



be such that P > 0. Thus, if k/e' > 0, p may have a limiting value of yk/e' and/or 
'k/e'. And if m 2 > e'(e 2 + g 2 ), other finite limiting values of p will also occur. If e' < 0, 



the range of p may extend to ±oo. 

One obvious advantage of the line element (|3U|) for the C-metric and its generalizations 
over more familiar forms is that it possesses well-behaved limits when a vanishes. 



22 



7 Non-expanding solutions 

The above sections have explored all the various cases of the Plebahski-Demiahski family 
in which the repeated principal null congruences have non-zero expansion. This has been 
achieved after first extending the metric (jSj) to the form (JT!?j) by including a shift in p. 
Our purpose now is to examine the alternative situation in which these congruences are 
non-expanding so that they are not covered by the metric (jHJ. And, since solutions with a 
non-expanding shear-free null geodesic congruence must be twist-free, the twist must also 
vanish in this case. The solutions obtained are therefore necessarily the type D solutions 
of Kundt's class. (An associated family of conformally flat solutions will be considered 
later in section 18.2)1 . 

Since these solutions are non-twisting, u cannot be a twist parameter and, in general, 
we may put u = 1. However, a special case occurs when u = 0. This will be treated 
separately in section 17731 below. 



7.1 Another modified form of the metric 



For the case in which w ^ 0, we start with the metric (|5jl . put u = 1, and apply the 
transformation 

r = 7 + nq, a = k~ x t, t = if> — ^ 2 kT x t. (43) 

where 7 and k are arbitrary parameters. (There is no need to include a shift in p at this 
stage. The possibility of including such a transformation will be considered later.) The 
resulting line element is then 



%{k^ - (l 2 +p 2 )dty ~^{dif) + (2 iq + Kq 2 )dtJ - P - dp 2 - £ dq 



P 



Pi 
P 



Q 



where 



Q = 1 — 0(7 + nq)p , 
p 2 = (7 + nq) 2 +p 2 , 

P = k + 2np - ep 2 + 2amp 3 - [a 2 {k + e 2 + g 2 ) + A/3]p 4 , 

Q = e + e x q - e 2 q 2 - (2cm + 4a 2 k>y + ^A<y)nq 3 - (a 2 k + |A)/t 2 g 4 , 



and 



K~ 2 (k + e 2 + g 2 - 2m7 + e-y 2 - 2cm7 3 - a 2 A;7 4 - l:A7 4 ) 



(44) 



e\ = 2KT x (-m + 67 - 3cm7 2 - 2a 2 fc7 3 - IA7 3 ) , 

e 2 = — e + 6cm7 + Qa 2 k^ 2 + 2A7 2 . 

We may make the simple observation at this point that the metric (jBj) explicitly includes 
both expanding type D solutions when k^O and non-expanding solutions when k — 0. 



7.2 Non-expanding cases with u 7^ 

To obtain the general family of non-expanding solutions, we now put k = in the above 
metric. (It can be seen from (JBJ) that, since Q{^) ~ K 2 eo here, the expansion and twist of 
the repeated null congruence must vanish in this limit.) For space-times to exist in this 
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limit, it is obviously required that eo and ei must remain finite. However, it is always 
possible to choose 7 to set t\ = 0. The metric is then 



d. 2 72 + p2 



Qdt 2 -l dq 2 - (y l p2)2 (ty + 2 1Q dt) 2 - i dp 



(45) 



(1 — cryp) 2 

where Q = e — e 2 g 2 and 

P = k + 2np - ep 2 + 2amp 3 - [a 2 (k + e 2 + g 2 ) + A/3]p 4 . 

It can be seen that the first terms in ()45|) represent a (timelike) surface of constant 
curvature whose sign is that of e 2 . It is then possible to use a linear transformation in q 
(and ip) together with a rescaling to set Q to one of three distinct canonical forms 

Q ■ l-g 2 , g 2 -l, 1. 

To achieve this, we require that 

e = -e 2 + 60717 + 2(3a 2 k + A)rf . 

In addition, 7 has been chosen to satisfy the equation 

3m + (e 2 — 2e)7 + 3cm7 2 = , 

This equation relates 7 closely to the parameter m. It is therefore appropriate to generally 
make the relabeling 7 = fh. We also have the constraint that 

k + e 2 + g 2 — mm + |(e + e2)m 2 = n 2 eo , 

in the limit as k — > 0. This effectively determines fc. (In the alternative forms for Q with 
different values of eo and e\, the coordinate g only covers a restricted range of the complete 
space-time. To transform between the canonical types and these others, a transformation 
of ijj is also required, as in the similar situation described in detail in subsection 15.3. 21 ) 

This family of solutions is thus characterised by one discrete parameter e 2 = 1,-1,0 
and six continuous parameters m, n, e, g, A and a. Of course, the range of p must be 
chosen such that P > with conformal infinity now at p — (am) -1 . But it is not now 
possible to make a linear transformation of p to simplify the roots of P without modifying 
the form of the metric ([45)1 . 

It can be seen that the conformal factor can be simplified by putting 

p — am 3 , , 

P = T— — . (46) 
1 + amp 

and applying a general rescaling. With this, P remains a quartic but with coefficients 
which are complicated combinations of the above parameters. 

This family of solutions was given by Carter j2] as his form [J3(— )]. It was also obtained 
by Plebahski who referred to these as anti-NUT metrics since the transformation (}4"3*|) 
is analogous to (|11|) which gives rise to generalised NUT metrics. For the vacuum case 
with A = 0, these solutions are of Kinnersley's class IV. An alternative form of the metric 
for this case was given by Kinnersley |4*Tj . 
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7.2.1 The case when a = 

For the case in which a = the metric (J45|) becomes 

ds 2 = p 2 (Qdt 2 -idg 2 ) - ^r(d^ + 2rhqdt) 2 - ^ dp 2 , (47) 
V Q ) P K ' P 

where p 2 = m 2 + p 2 , P = k + 2np — ep 2 — |Ap 4 in which e = 2 Am 2 — e 2 , and m is 
chosen to satisfy the equation |Am 3 — e 2 m — m = 0. It is also necessary that k is given 
by A; = — e 2 — g 2 — e 2 m 2 + Am 4 + eo/t 2 , which is taken in the limit as k — ► 0. 

If m ^ 0, these solutions are nonsingular. The function P is a quartic if A 7^ 
(otherwise it is at most a quadratic). It may therefore have up to four roots and the 
range of p must be chosen such that P > 0. The metric for these solutions is now given 
by ()47|) in which Q takes one of its canonical forms, p 2 = p 2 + m 2 , and 

P = -(e 2 + g 2 + e 2 m 2 - Am 4 ) + 2np + (e 2 - 2Am 2 )p 2 - §Ap 4 . (48) 

When the parameter m, the charge parameters e and g, and the cosmological constant 
all vanish (with n ^ 0), these solutions are the B- metrics [40J in different coordinates: 
BI when t 2 — 1? BII when e 2 = —1, and Bill when e 2 = 0. When these additional 
parameters are non-zero, they represent a large family of distinct space-times. These 
depend, not only on the canonical form of Q, but also on the possible roots of P and the 
range of p that is chosen to ensure that P(p) > 0. 

For the case in which e 2 = 1 (i.e. when Q = 1 — q 2 ), we can apply the transformation 

p = r, q = cos 9, ip — —(p. (49) 

This gives the line element 

ds 2 = (r 2 + m 2 )(sm 2 6dt 2 - dd 2 ) - — — f d(j) - 2m cos 9 dt] — dr 2 (50) 



where 



(e 2 + g 2 + m 2 - Am 4 ) + 2nr + (1 - 2km 2 )r 2 - |Ar 4 . 



When m = e = g = A = 0, n^0, this is the more familiar form of the B /-metric. 
However, it is not necessary or desirable to transform q so that it is restricted to the 
range q G [—1, 1]. It is better to retain a general form which can be extended through the 
horizons at q = 1 and q = — 1. 

Since the repeated principal null congruence is non-expanding and non-twisting (as 
well as being geodesic and shear-free), this family of solutions must belong to the type D 
solutions of Kundt's class (see section §31.7 of [IS]). Indeed, the metric (}4Tj) can be 
written in the standard form of the Kundt solutions using the transformation 



/ f 1 t f d<? 

p, y = ip + 2m J --^dq, u = t ~J~Q> 



P 2 Q, 



and putting 

V2( = x + iy, x = Jp 2 (z)dz, 



where 



P' 



This leads to the metric 



ds 2 = 2du(dv + W d( + W d( + H du) - 2d p 2 ^ , (51) 
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where 



z 2 + m 2 



W 



H 



-(e 2 + g 2 + e 2 m 2 - Am 4 ) + 2nz + (e 2 - 2Am 2 )z 2 
y/2v 



lAz 4 



(z + zm)P 2 ' 
\e (z 2 + m 2 ) - 



?2 



+ 



2m 



_2(z 2 + m 2 ) (^ 2 +m 2 ) 2 P 2 _ 

This is a generalization of the solution given as equations (31.41) and (31.58) in which 
now includes a cosmological constant and an additional discrete parameter. Further, there 
is a different identification of the parameters m and n. 



7.3 Non-expanding cases with uj = 

Let us now consider the remaining exceptional case in which we again start with the 
metric (0), but now with u = 0: i.e. we start with the metric ()39|) with ([40)1 . We then 
construct non-expanding solutions for this case by applying the transformation 

7 2 

r = 7 + nq, t = — t, 

K 

and taking the limit in which k = 0. This procedure leads to the line element 



ceyp)' 



Qdt 2 -ld q 2 -Pda 2 - 1 d p 2 
Q P 



where 



and 



P = k + 2n'p — ep 2 + 2amp 3 — a 2 (e 2 + g 2 )p 4 , 
Q = e + e iq - e 2 q 2 , 

e = fiT 2 (e 2 + g 2 — 2m7 + e7 2 — 2an'7 3 — a 2 A;7 4 — |A7 4 ) , 
Cl = 2KT x {-m + €7 - 3cm'7 2 - 2a 2 ky i - f A7 3 ) , 
e 2 = — e + 6an'7 + 6a 2 /ry 2 + 2A7 2 . 

Thus, Q is again an arbitrary quadratic which, for finite coefficients determines the ap- 
propriate choice of 7 and provides a constraint on the other parameters. 

The case in which a = will be dealt with in the following section. However, when 
the conformal factor can be removed by the transformation 

1 — a^p = ^ pT x . a = aa. 



This takes the metric to the form 

ds 2 = p 2 (q dt 2 - 4 dg 2 ) - ^-da 2 - 4 dp 2 , (52) 



Q J P P 



where Q may take one of the three standard canonical forms and P(p) is a quartic function 
whose coefficients are different combinations of the parameters m, n', e, k, e, g and 7 
to that given previously. This metric is in fact equivalent to (|45|) after applying the 
transformation ()46|) in the singular case in which 7 = (i.e. m = 0) in the notation of 
that section. However, the coefficients of the quartic P are different. 
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This family of solutions has a curvature singularity at p = 0. However, since P is 
generally a quartic and we require that P > 0, different ranges of p are possible. The 
different possibilities which arise therefore represent various special cases of both singular 
and non-singular space-times. 



8 Non-twisting and non-accelerating solutions 

8.1 The metric (03) with a = uj = 

Let us finally consider the case when both the acceleration and the twist vanish. Putting 
a = uj = accordingly, the metric © becomes 

ds* = %di* - ^dr 2 - t* (±-dp> +Pd<A , (53) 
r 2 Q \P J 

where 

P = k + 2n'p - ep 2 , Q = (e 2 + g 2 ) - 2mr + er 2 - ±Ar 4 , (54) 

as can also be seen from (J3l?j) and (|4Up. Apart from the possible presence of a non-zero 
cosmological constant and hence curvature scalar, the only non-zero components of the 
curvature tensor take the form 

m e 2 + g 2 ^ e 2 + g 2 

^ = --3+—^' $11 = ^^' (55) 

since the parameter n must vanish in this case. 

The surfaces on which r and r are constant have constant curvature e, so P may be 
set to one of the standard three canonical forms. For the case in which e = 1, P = 1 — p 2 , 
and we can put 

p = cos9, a = —(f), T — t, 

and the metric (}53|) then takes the familiar form of the Reissner-Nordstrom-de Sitter 
solution 



, , , 2m e 2 + g 2 A 2 \ 2 / 2m e 2 + g 2 A 2 \ 1 , 2 

ds 2 =(l + T~ + -r 2 dt 2 - 1 + ^f— + — r ) dr z 

r t 3 / \ r r z 3 

-r 2 (d^ 2 + sin 2 0( 



This obviously includes the Schwarzschild, Reissner-Nordstrom, de Sitter and anti-de Sit- 
ter solutions and their various combinations. Of course, this is the particular nontwisting 
subcase of the metric discussed in section ll~2l (for a = = I). 

Associated, but less physically significant, cases occur when e = (P = 1) and when 
e = — 1 (P = p 2 — 1). For the vacuum case with A = 0, the different permitted values of 
e cover the three versions of the so-called A-metrics 



8.2 Extension to include non-expanding solutions 

The above family of solutions covers all non-accelerating and non-twisting type D solutions 
with non-zero expansion. Interestingly, it is also possible to perform a coordinate trans- 
formation of which, in a certain limit, includes an associated family of non-expanding 
type D and conformally flat solutions. 
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Starting with the metric ()53|) with (|54p. let us consider the coordinate transformation 
which involves shifts in both r and p simultaneously 

p = (3 + Kp, r = 7 + Kg, a = k~ a, t = 6 2 /t _1 f, (56) 

where /3, 7, k and b are arbitrary constants. With this, the metric (J53|) becomes 

d * 2 = (7^|)i df 2 " (i + ^ 2 (g d ? + ? * 2 + ™ s ' 2 ) • < 57 > 

where 

P = a + aip + a 2 p 2 , 

Q = 6 + hq + 6 2 g 2 - ^Al^q 3 - ±A/t 2 g 4 , 

and 

62 = e - 2A7 2 a 2 = -e 

61 = 2«r 1 (-m + C7 - |A7 3 ) 0l = 2ftT 1 (n' - e/3) 

6 = K " 2 (e 2 + g 2 - 2m 7 + e 7 2 - |A 7 4 ) a = /€" 2 (A; + 2/5n' - e/3 2 ) 

As in the previous section, we proceed to the non-expanding case by setting k = 0. In 
this case, the metric becomes 

ds 2 = b 2 (Y df 2 - i dg 2 ) - 7 2 (P da 2 + i dp 2 ) , (58) 

where Y = 6 2 7 _2 (6o + b\q + ^g 2 )- Provided 62 and e are non-zero, the linear terms in Y 
and P can be removed by choosing the free parameters 7 and j3 respectively such that 

|A 7 3 - e7 + m = and e/3 - n' = 0. (59) 

In addition, the constant terms in Y and P will only be bounded if the parameters of the 
solution satisfy the constraints 

e 2 + g 2 - 2my + e7 2 - iA7 4 = b K 2 . 

3 (60) 

k + 2(3n' - e(3 2 = a n 2 , 

as n — > with b and a bounded. Together with (JHHJ), the second of these constraints 
determines the required value of k for such solutions to occur, and the first constraint is 
a condition on the parameters m, e, g, e and A. 

In all cases, linear transformations in p and q can then be used, together with a 
rescaling, so that P and Y respectively each take one of the standard canonical forms. In 
this way, we obtain the family of solutions given by Carter jH] as his form [D] , which has 
also been given by Plebahski as his case C. 

From the above discussion, together with that given below, it can be seen that the 
metric (|57j) explicitly contains both the family of Schwarzschild-Reissner-Nordstrom- 
de Sitter space-times (when k ^ 0) and the space-times that are a product of two 2- 
dimensional spaces of constant curvature such as Bertotti-Robinson (when k — 0). For 
the case in which the cosmological constant vanishes, the solution (|57j) includes that of 
Ray and Wei [121 (rediscovered by Halilsoy jSj) which combines the Schwarzschild, the 
Reissner-Nordstrom and the Bertotti-Robinson solutions in a single metric. 
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8.2.1 The Bertotti— Robinson solution 

Let us now investigate the metric (J58j) for the case in which A = and g = 0. In particular, 
for the case in which e > 0, we use a rescaling to put e = 1. We may also choose 6 = 7 
and then choose j3 = n' and 7 = m to satisfy (J5T?|) . We then consider the case in which 
the parameters satisfy the constraints ()60|) with k = —n' 2 + k 2 in the limit as k — > 0. 
In this case P = 1 — p 2 and Y = q 2 . In addition, we have e 2 = m 2 so that we are 
now considering a degenerate limit of the extreme Reissner-Nordstrom solution. In this 
case, the Weyl tensor component \&2 vanishes identically. The resulting space-time is a 
conformally flat non-null electrovacuum solution, which must be the Bertotti-Robinson 
solution. Indeed, with the coordinate substitutions p = cos 6, a = —(f) and q — 1/f, the 
line element for this case becomes 

e 2 

ds 2 = — (df 2 - dr 2 - r 2 d6 2 - r 2 sin 2 9 d0 2 ), (61) 
which is the most familiar form of the Bertotti-Robinson solution. 

8.2.2 Other direct product space-times 

For the remaining cases that are included in the metric ()58|) with the conditions ()59|) and 
(|60|h we note that P can always be put into one of its three distinct canonical forms 
which correspond to the different parametrizations of 2-spaces of constant curvature. In 
addition, with 7 chosen to satisfy (}59|) . Y is given by 

Y = e - e 2 q 2 , 

where e 2 is given by 

e 7 ~ 2 + e 2 r 2 = 2A, 
and eo satisfies the constraint ()6U|) which becomes 

7~ 2 (e 2 + # 2 ) -e + A 7 2 = K 2 e b' 2 = 0. (62) 

In fact, it is possible to choose b so that e 2 is equal to +1, —1 or 0. Then, since 7 is already 
defined by (J5T?|) in terms of m, e and A, (Jo^j) represents a constraint on the parameters 
that is required to avoid a divergence in the metric. 

The non-zero components of the curvature tensor (apart from the Ricci scalar which 
is proportional to the cosmological constant) are now 

* 2 = -|A, $ u =I 7 -V + <7 2 ). 

(The terms containing m and e 2 + g 2 that occurred in ty 2 in f|55|) have cancelled out in 
the limit n = 0, apart from a term proportional to the cosmological constant, through 
the choice and constraint (j5^|) and (JH21)-) 

It can thus be seen that the metric (|58|) is the direct product of two 2-dimensional 
spaces of constant curvatures e and e 2 of signatures (+, — ) and (— , — ) respectively. The 
values of these two parameters determine entirely the geometry of the space-time. As 
described e.g. by Ortaggio and Podolsky [H], these represent the complete family of 
Bertotti-Robinson, Narai and Plebahski-Hacyan space-times. Excluding cases with neg- 
ative energy density, it is found that there are six possible geometries. 

Allowing for the possibilities of both electric and magnetic charges, this family of 
space-times essentially depends on only three free dynamical parameters, which may be 
taken as e, g and A. (The parameter m, which is necessarily non-zero, may be reintroduced 
through the definition of 7 and the constraint (JH21)-) 
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9 Conclusions 

We have analysed in detail the complete family of Plebahski-Demiahski metrics and 
those that can be derived from them by taking limits after performing some specific 
coordinate transformations. These non-null solutions of the Einstein-Maxwell equations 
with a possible cosmological constant include all type D solutions in which the principal 
null directions of the electromagnetic field (if included) are aligned with the repeated 
principal null directions of the Weyl tensor and the group orbits of the Killing vectors are 
non-null. 
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Figure 2: Vacuum type D solutions with A = and m ^ 0. In the expanding case, these 
solutions generally have four continuous parameters m, I, a and a, and two auxiliary parameters 
e and k that can be scaled to any convenient specific values. The parameters associated with each 
solution are shown above each box - those in brackets have discrete values. The abbreviation 
"gen" indicates that the named solution is generalised to include the associated solutions with 
other values of the discrete parameters. The Kinnersley class of each sub-family is also indicated 
below each box. All particular case have obvious generalizations which include charges and/or 
a cosmological constant. 



In general, we have found it convenient to treat separately the cases in which the 
repeated principal null congruences are either expanding or non-expanding. For the ex- 
panding case, the metric which includes all possibilities is given by (|12p. For the non- 
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expanding case, the general metric is given by although, when A ^ 0, other type D 
solutions are given by (JESJ), which also includes certain associated conformally flat solu- 
tions. The non-expanding solutions can alternatively be expressed in the standard form 
for solutions of Kundt's class, namely (f5*Tj). 

The structure of the complete family of Plebahski-Demiahski solutions for the vacuum 
case with e = g = A = 0is illustrated in figure El For this case, it seems to be most 
convenient to represent the expanding solutions in terms of four continuous parameters 
m, I, a and a, and two auxiliary parameters e and k that can be scaled to any convenient 
specific values. For the case in which the line element has a component in which the 
coordinates span surfaces of positive curvature, the classification of these solutions has 
been illustrated in more detail in figure ^ In that case, the family of Robinson-Trautman 
type D solutions is represented by the C-metric. It was appropriate there to distinguish 
the cases in which the Kerr-like rotation parameter is greater or less than the NUT 
parameter, as these cases have significantly different singularity structures. However, this 
distinction is not continued in figureEl Here, the solutions are generalized to include those 
with surfaces of alternative curvatures and the non-expanding solutions are also included. 
All the particular cases illustrated in figure El have obvious generalizations with non-zero 
charge parameters and/or a non-zero cosmo logical constant. 

In general, these solutions are characterized by two generally quartic functions whose 
coefficients are related to the physical parameters of the space-time. However, although 
it is traditional to use an available coordinate freedom to remove the linear term in one of 
these quartics, it has been found to be much more helpful to use this freedom to simplify 
the roots of at least one of these functions. This significantly simplifies the calculations 
involved in interpreting these solutions and in identifying the physical interpretation of the 
coefficients. In particular, as also shown in [21], it is found that the Plebahski-Demiahski 
parameter n is not the NUT parameter, but is related to it by (fT5|). 

Of course, most of the solutions in this family do not yet have any known physical 
significance as models of the gravitational fields of realistic sources. Nevertheless, they 
provide a most important family of model space-times that can at least be interpreted 
geometrically. 
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